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Abstract 

The aim of this paper is to show, that the 'oscillating universe' is a viable alternative 
to inflation. We remind that this model provides a natural solution to the flatness or 
entropy and to the horizon problem of standard cosmology. We study the evolution 
of density perturbations and determine the power spectrum in a closed universe. The 
results lead to constraints of how a previous cycle might have looked like. We argue 
that most of the radiation entropy of the present universe may have originated from 
gravitational entropy produced in a previous cycle. 

We show that measurements of the power spectrum on very large scales could in 
principle decide whether our universe is closed, flat or open. 



1 Introduction 



In a closed universe, the inevitable big crunch might actually be followed by a subsequent big 
bang. This idea is very old. It goes back to Lemaitres Phoenix picture of 1933 As we shall 
remind especially younger readers, this 'oscillating universe' provides quite naturally a solution 
to the flatness or entropy and to the horizon problem. This was well known before the advent 
of inflationary models around 1980 in seminal papers by Starobinski, Guth, Linde and others 
0. Since then, it has been so completely forgotten, that we ourselves originally believed to 
be the first to study these ideas and older colleagues had to refer us to the original literature. 
However, it is not only inflation which let people forget about the oscillatory universe, but 
also an argument in essence due to Penrose that black hole formation in a previous cycle 
leads to far too much entropy and thereby to an even more severe entropy problem in the 
opposite sense than the usual one. 

It is the aim of this paper to show that under most realistic circumstances Penrose's 
conclusion need not be drawn. We shall see that some amount of entropy production due to 
gravitational clumping can just somewhat accelerate the growth of the maximal scale factor, 
c^max, from one cycle to the next, without over producing entropy. This will lead us to the 
conclusion that the oscillating universe remains a viable alternative to inflation. 

We consider this especially important since inflation has become some kind of 'cosmological 
dogma' during the last ten years, despite the fact that no inflationary scenario which solves 
the horizon and flatness problems and yields acceptably small density fluctuations has yet 
been constructed without substantial fine tuning (which may be protected by a symmetry 
and thus be 'technically natural'). Furthermore, many inflationary models are built upon the 
gravitational action of vacuum energy, which acts gravitationally like a cosmological constant, 
the most miraculous number in cosmology, which today is by a factor of about 10^°° times 
smaller than what we would expect from particle physics 0]. This unbelievable amount of 
fine tuning tells us, that our understanding of the interplay of gravity and quantum vacuum 
energy is unsufficient. Therefore, a mechanism relying mainly on this interplay, to us, seems 
unsatisfactory. 

Besides solving the horizon and entropy problems, inflation generically predicts a scale 
invariant Harrison-Zel'dovich initial spectrum of fluctuations as it was observed by the DMR 
experiment on the COBE satellite These observations have been considered as great 
success of even 'proof of inflation. However, also global topological defects or cosmic 
strings which can form during phase transitions in the early universe, naturally lead to a 
scale invariant spectrum of fluctuations but they cannot easily be reconciled with inflation. 

These considerations prompted us to look for possibilities to solve the flatness and the 
horizon problem without invoking an inflationary period. 

The basic picture which we work out in this paper is the following: 
The first 'big bang', the first 3-dimensional closed universe, emerged from quantum fluctuation 
of some, e.g., string vacuum. Its duration was of the order of a Planck time. Due to some 
non thermal processes there may have been a small gain of entropy, S^^ < S^lj^. The first 
big crunch triggered the formation of the next big bang who's entropy was slightly larger and 
therefore its duration was slightly longer, S^^ < S^^^^ < S'^ . This process continues with ever 
longer cycles. In a cycle with duration significantly longer than Planck time, we assume, that 
after a few Planck times during which the universe may have been in some quantum gravity 
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or stringy state, we have a mainly classical, radiation dominated universe. With the exception 
of short periods during which matter and radiation fall briefly out of thermal equilibrium, the 
universe expands and contracts adiabatically. 

Since we have no theory for the string or quantum gravity era of cosmology, it is very 
difficult to estimate the number of Planckian cycles before the universe enters a radiation 
dominated era. Furthermore, the entropy producing processes at high temperature depend on 
the microphysics at very high energy. We thus do not have the tools to estimate the number of 
cycles which never entered a low temperature era. But, as we shall see, due to the production 
of gravitational entropy, the number of cycles which entered a matter dominated era is very 
limited. 

As long as the universe remains radiation dominated no black holes can form. Only during 
a cycle with a long enough matter dominated era black hole formation is possible. Penrose's 
argument now goes as follows: During a matter dominated era (small) black holes can form. 
These finally, during the collapse phase, coalesce into one huge black hole which eventually 
contains the whole mass of the universe. Its entropy is thus given by 

= (1/2) = 2nRl/G = 8nGM^, > 10'^\ (1) 

where we have set Mbh equal to the present mass within one Hubble volume which is of order 
10^^ Mq to obtain the last inequality. Clearly, already a significantly smaller mass would 
do, since the actual radiation entropy within the present Hubble radius is of the order of 
S Hubble ~ 10^'', a discrepancy of nearly 40 orders of magnitude. In terms of entropy per 
baryon, this yields r]~^ = S/Nb ~ 10"^^ instead of the observed value ri~^ ~ 10^. 

Is there a way out of this simple but disastrous conclusion? 
The first objection is that the laws of black hole thermodynamics which rely heavily on 
Hawking radiation, hold only in asymptotically fiat spacetimes. Or, at least, that the entropy 
of a black hole can be set equal to its area only for an observer outside the black hole itself. 
Therefore, the black hole entropy formula should only be adopted for black holes much smaller 
than the size of the universe. Let us therefore add only the entropy of black holes which are at 
least 10 times smaller than the curvature radius of the universe and neglect subsequent growth 
of entropy due to the coalescence of these black holes. Of course, this rule is somewhat ad hoc, 
but as long as we have no clue of how to calculate in general the entropy of the gravitational 
field it is a possible 'rule of thumb'. However, with this correction we gain only about a factor 
of 10 in the above entropy formula (|l]) and not the required factor of about 10^^. But there is 
an additional short come in the Penrose conclusion: The radiation entropy which we observe 
today is the entropy generated mainly during the previous cycle whose matter dominated 
epoch might have been much shorter, leading to significantly less clumping and thus much 
less gravitational entropy production. 

Furthermore, as we shall see, it is not clear that structure forms via hierarchical clustering. 
In a pure radiation universe, it may well be that large black holes form first (if at all!) and 
the black hole entropy formula cannot be applied. 

From these arguments it should be clear, that Penrose's objection to the oscillating universe 
does not have to be accepted and there may be ways out. Another possibility not investigated 
in this work is Israel's idea of mass infiation inside the horizon of black holes p], ^. Israel et 
al. accept the black hole entropy formula, but argue that inside the black hole horizon mass 
infiation takes place such that the ratio rj^^ gets reduced substantially. 
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The reminder of this paper is organized as follows: In the next section we give a brief 
review of the oscillatory universe. In the main part of this paper, section 3, we investigate 
cosmological perturbation theory in a closed universe and determine the evolution of a Harrison 
Zel'dovich initial spectrum in a purely radiation dominated universe and in a universe with 
an intermediate matter dominated epoch. We discuss how the final spectrum depends on the 
duration of the matter dominated epoch and we formulate a limit for the maximum radiation 
entropy of the previous cycle. The final section is devoted to our conclusions. 

Notation: The scale factor of the Friedmann universe is denoted by a, we use the con- 
formal time coordinate and the metric signature (— , +, +, +), so that the Friedmann metric 
is given by 

ds"^ = a^{dt^ - ^ijdx'dx^) , (2) 
where is the metric of the unit three sphere, e.g., 

'jijdx'dx^ = dx^ + sin^(x) (c?6'^ + sin^ ddip"^) . 
Cosmic time is denoted by r, r = /* adt. 

We normalize the scale factor a such that the curvature of the spatial sections is equal to 



2 Reminder to the oscillatory universe 

Let us first explain how the oscillatory universe solves the flatness or entropy problem. To 
do this it is useful to state the problem in a somewhat different form: If the equation of 
state in a Friedmann universe satisfies the strong energy condition, p + 3p > 0, then p decays 
faster than and = 1 is the unstable 'initial' fix point of expansion. This means each 
Friedmann universe starts out at t ~ a few ^pja^ck ^i^^ i7 ^ 1 and later deviates more and 
more from this value. The flatness problem can thus be stated as follows: How can it be, 
that our universe at its old age, t ^ ^pianck' ^ -^Planck ^^^^^ looks so young, ~ 1? This 
problem is easily solved in the oscillating universe as we shall now show. 

The following arguments are due to Tolman [l^ . Only a year after Lemaitre first brought 



up his phoenix picture Tolman realized: Since the entropy of the next universal expansion can 
only be larger than the previous one, the maximum expansion factor of the next cycle, amax, 
is larger than the corresponding maximum of the previous cycle. Since the density parameter 
Q starts deviating from 1 only when the scale factor a approaches Omax, in the next cycle it 
will take longer until this happens. We consider now a cycle with a duration substantially 
longer than Planck time which has entered a radiation dominated phase. If relativistic matter 
is in thermal equilibrium (which we assume to be true most of the time) its energy density 
and entropy density are given by = c = ^Boltzmann ~ ■'-) 

p = —NT'' (3) 
30 ^ ^ 

s = ^NT', (4) 
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where denotes the effective number of degrees of freedom (spin states). N = Nb + (7/8) A^/. 
Here Nf, are bosonic degrees of freedom and Nf are fermionic degrees of freedom. Furthermore, 
from Friedmann's equation for a closed universe, 

+ 1 = P , (5) 



,a/ 3 
together with (|]) and (0) one finds 

p- pc SnGpa'^ - 3(a/a)2 1 

i2 — 1 = = = -j:, . [b) 

Pc 3{a/ay Q (^i^y' S'^/^T'^ - 1 

Here 5* is the total entropy of the universe, S = 27r^a^s. Therefore, the larger the total 
entropy S the smaller the deviation of Q from the critical value 1 at a given temperature T, or 
the lower temperatures are required for a substantial increase of Q. Expressing the maximal 
scale factor, Omax, the age of the universe at maximal expansion, t^^x = T(amax) and the 
minimal temperature, T^in = 7'(amax), in terms of the entropy also show that these values 
grow, respectively decrease with increasing entropy: 

adt = Omax , imax = 7r/2 , (8) 

r„. ^ A,s-v^ , A,^(j-|^)"°. (9) 

The time it takes for the density parameter to differ significantly from 1 is a substantial 
fraction of Tmax- Therefore, the universe 'looks young' for longer and longer times as the 
entropy increases cycle by cycle. 

It is clear that in the oscillating universe also the horizon problem disappears since the 
age of the universe is not given approximately by the inverse Hubble time, which is the age 
of the present cycle, but the sum of the ages of all previous cycles has to be added, leading to 
a much larger age which might even be infinite. For this solution to be valid, it is important 
that correlations are not lost during a big crunch/big bang passage and that the behavior of 
particles or strings during this time is governed by a causal theory. A similar problem is also 



encountered in 'pre-big-bang' string cosmology [[11 1. 

In Appendix B, we explore the possibility that quantum gravity may effectively lead to an 
Euclidean region of spacetime close to the big crunch/big bang era. This example of a causal 
continuation from one cycle to the next is due to Ellis [0, |1^. The singularity in the metric 
induced by the signature change is very mild. We show how in this case geodesies can be 
continued through the crunch in a completely smooth manner. 



3 Cosmological Perturbation Theory in a Closed Uni- 
verse 

In this section we first study the equations which govern the time evolution of radiation and 
matter density perturbations in a closed universe. We then determine the power spectrum. 
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which we use to decide at which length scale the density perturbations might first lead to the 
formation of objects (e.g. galaxies or black holes). We finally use these results to argue how 
the entropy of the present cycle may have been generated. 



3.1 Time Evolution of Density Perturbations 

To describe the time evolution of density perturbations we use gauge invariant linear cosmo- 
logical perturbation theory (see e.g. [0). Assuming adiabatic perturbations and neglecting 
anisotropic stresses, the evolution of the gauge invariant density perturbation variable D is 
governed by the equation 

i) - (V^ + 3)c2D + (1 + 3c2 - Qu) b (10) 

-3 1^ - 3 ^""^ - ^) + (1 + ^)\'^Gpa'^ D = 0. 

In a universe which consists of matter and radiation, u = p/p = {l/3){l+a/aeq)~^, = p/p = 
(1/3)(1 + 3a/4aeg)~^, where Cg is the sound velocity. A dot indicates derivatives with respect 
to conformal time t and Qeq is the scale factor when prad = Pmat- Two cases of particular 
interest are dust (t^ = = 0, = 0) and radiation (t^ = = 1/3, a^q = oo). 

Expanding D in terms of scalar harmonic functions on S'^, as described in Appendix |^, 
leads to the following equation for the gauge invariant density perturbation amplitude for the 
wavenumber / G {0, 1, 2, ...}: 

biit) + (/(/ + 2) - 3)clDi{t) + (1 + 3cl - Qu) biit) 

3 {^^\cl + + u)^^Gpa'\ D,{t) = 0. (11) 



a 

UJ I - 

a 



For most of the sequel we omit the index / which distinguishes the different eigenfunctions of 
on S^. 

In the following subsections we solve equation (|l^) in some cases of special interest. We 
then use our results to derive the power spectrum. 

3.1.1 Radiation Density Fluctuations in a Radiation Universe 

At early stages of expansion and at the end of the collapsing phase, the universe will consist 
of pure radiation, i.e., all matter will be relativistic. Therefore, this case is important for each 
hypothetical previous cycle, whether it entered the matter dominated era or it was always 
radiation dominated. For radiation, where = = 1/3, equation ( |Tll) reduces to 

'' (/(/ + 2) -3) fd\ 16 ^ 2^ 



D + <i \ -[-)- yTTGpa' \D = Q. (12) 
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Figure 1: The time evolution of Di{t) for radiation density perturbations in a radiation 
dominated universe (in arbitrary units). The amphtude Di{t) is shown for the scales / = 20 
and / = 80. Very close to the big crunch, the divergence due to cot t takes over. 



Inserting the solution of the Friedmann equation for the scale factor of a radiation dominated 
universe, which is a{t) = OmaxSint with amax = (SvrGpa'^/S)^/^, we obtain 



(sin^ t)D + 



'/(/ + 2) 



sm^t-2]D = 0, te[0,7i]. 



(13) 



The solution of this equation is given by 



D{t) = sinh - 



1 d 



sin t dt 



[ci exp{i^/a'lt} + C2 exp{-iy/ait}] 



(14) 



with ai = l{l + 2)/3 and / ^ (see [0). We are only interested in real solutions for positive 
integers Z, in which case D{t) can be written in the form 



ci [^/al cot t ?,m.{^/a'lt) — ai cos{^/alt)] 
+ C2[^/(^l cott cos{^/alt) — aism{y/ait)] . 



(15) 



This solution is plotted in Fig. |I]for / = 20 and / = 80. Obviously the amplitude cott diverges 
at the big bang and at the big crunch where t = and t = n respectively. Since we assume 
that the fluctuations are created at some time tj > after the big bang, the divergence at 
t = is not a problem. Apart from its oscillation with frequency ~ /, the amplitude of 
density fluctuations is approximately constant for most of the cycle, but diverges close to the 
crunch like Di(t) oc /(vr — t)~^. 
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3.1.2 Matter Density Fluctuations in a Radiation Dominated Universe 

For dust (u; = Cg = 0) equation (0) reduces to 

D + -D- AirGp-raata^D = 0. (16) 
a 

As long as the universe expands, d is positive and hence the second term in this equation acts as 
a damping term. This term vanishes at maximum expansion and turns into a stimulation when 
the universe contracts. Therefore we expect the growth of fluctuations to become substantially 
enhanced during the contraction phase. 

Inserting the scale factor a{t) = amax sin t for the radiation dominated universe, we obtain 

sintiD + costD - fiD = 0, te[0,7r], (17) 

where fx = AnGprnatd'^ = (3/2)(amaa:/aeg) and Umax = (STrGpa^/S)^/^. With the substitution 
a; = sint this equation leads to 

x{x - l)(a; + 1)D" + {2x^ - 1)D' + fiD = 0. (18) 

This equation is a special case of Heun's differential equation. For < 1 one solution 

oo 

Daix) = l + Y^ CnX"" 
n=l 

of (|T8|) is a convergent power series with 

Ci = /i 

_ n(n-l) 

" (71 + 1)2''"-'+ (n + 1)2 



A numerical solution of ([T7| ) is shown in Fig. |^. 

Let us discuss the behavior of matter density fluctuations more closely during the different 
epochs of a cycle. At early times, when t <^ 1, equation (plTl) simplifles to 

W + D - fiD = (19) 

which has a solution in terms of Bessel functions: 

D{t) = ci Jo(2zv7^) + C2Yo{2iy/JIt) ^ di + C2ln{2y/Jd) . 



Neglecting logarithmic growth, these fluctuations are approximately constant for t ^ 1. Fig. ^ 
shows, that the logarithmic growth for small t is a good approximation up to t ~ 1/10 which 
means that all scales with / > 10 or so enter the horizon during this era, where we can consider 
the fluctuations to be approximately constant. This is the well known Mezaros effect: matter 
fluctuations do not grow in a flat radiation dominated universe. Close to maximum expansion 
(t = n/2 ± e, |e| ^ 1), equation ( p!7|) reduces to 

d^D dD 

-^-e fiD = 
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Figure 2: The time evolution of D{t) for matter density perturbations in a radiation domi- 
nated universe (in arbitrary units). For comparison, a logarithmically divergent fit is shown 
close to the coUaps (solid line). The vertical dotted hne indicates t^ax, the time when the 
cycle reaches its maximum expansion. 



The solution can be written in terms of the confluent hyper geometric function: 

/U 1 1 —e 
2 ~ 4'4'^ 



Die) = e-^V/^Y(^-ii^] , (20) 



with 

Y{k,m,x) = CiMk^m{x) + C2MA. _„(x) 

Mfc,„(x) = x^/^+'^e-^/2F(l/2 + m-k,2m + l,x) 

Near maximum expansion, small |e|, D grows nearly exponentially, like in a non-expanding 
universe. 

Close to the big crunch, where t = tt — e, we again obtain equation ([19|) , replacing t 
and dots by e and derivatives with respect to e. This reflects the symmetry of the closed 
universe between (big bang, t) and (big crunch, —t) as long as the entropy remains unchanged. 
Therefore, close to the big crunch D diverges logarithmically: 



D{t) oc In 



1 



This solution is valid until the particles become relativistic around the de-confinement phase 
transition where T = 100 MeV, from where on we have to consider pure radiation density 
fluctuations. 
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3.1.3 Matter Density Fluctuations in a Matter Dominated Universe 

In the case of a matter dominated universe, the scale factor is given by a{t) = {a^ 
cost), where a^ax = SuGpa^/S and hence (p!6| ) reads 



(1 - cost)D + (sint)Z) - 3D = 0, 
with the well known solution (see |]T6| ) 
5 + cos t 3t sin t 



D{t) = ci 



cost 



cost)^ 



+ C2 



sint 



cost)" 



(21) 
(22) 



At early times, equation ( PT| ) is approximately given by t'^D + 2tD - 6D = 0, such that 
D cct^ or D oc t~^, for the growing and decaying mode respectively. Close to the collapse, for 
t = TT — e, equation (pT]) reduces again to 



e^-^D + 2e—D-QD = 0, 
cte^ ae 

but the growing and decaying modes are interchanged. Now the growing mode solution is 
given hj D = Doe~^ = Dq{2it — t)~^. At maximum expansion, the damping term vanishes 
and the evolution of the fluctuations around t = tt is described by exponential growth or 
decay. 



3.1.4 Composite Model 

To construct a more realistic model where the scale factor is not only determined by a single 
matter or radiation background, we now assume a simple composite model, where the energy 
density of the universe is given by 



p(a) 



Peg 

2 



"eq 



+ 



"eq 



The first expression on the right hand side represents the behavior of the matter den- 
sity and the second term reflects the behavior of radiation density. The solution of the 
Friedmann equation (|^) in this case is 



a(t) = sin ( t — arcsin 



+ a = a^qia sin t + sin^ -) 



(23) 



where a = fleg/oo with ao = {4:7iGpeq/S) and 



^a'^tteq, A = a^a^g + a 



.,2^2 



7^2 



Furthermore, we find from (E3D that 



arcsin(A 



-1/2/ 



"eq 



a)) + arcsin(aA 



-l/2x 



TT 



/2 + arcsin(aA"^/^), 



a{tr 



(a + y/ITc^) 
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Figure 3: The amplitude D{t) of matter density perturbations in the composite model 
(in arbitrary units). The solid line shows D{t) for a cycle with a short matter dominated 
phase (a = 4). The dashed line shows D{t) for a cycle with a long matter dominated epoch 
{a = 1000, ttmax ^ CLeq)- The left and right vertical dotted lines indicate the time of maximum 
expansion of the cycle for a = 4 and a = 1000 respectively. 

We can use a as a parameter which determines the duration of the matter dominated epoch 
in a closed universe containing matter and radiation. For a -C 1, Omax ~ CiCLeq ^ CLeq and the 
universe never becomes matter dominated. For a 3> 1, Omax ~ «^fleg ^ CLeq and t^g <^ tmax ~ 
vr; the universe experiences a long matter dominated era. 

For radiation density perturbations in this composite model, equation ( |Tl| ) yields 

a''{t)b{t) + |a2(t) |^ ^(^ + ^) _ ^(^)^(^) _ 4^^^^! ^(^) = g, (24) 

and for matter density perturbations we obtain 

a{t)D{t) + d{t)D{t) - 3aD{t) = 0, (25) 

with a{t) given by (^5]). We are particularly interested in equation (P3|). Numerical solutions 
for the growing mode of a short and long matter dominated phase are shown in Fig. |^. 

3.2 The Power Spectrum 

3.2.1 The Harrison ZeI'dovich initial spectrum 

The power spectrum P{l,t) determines the scaling behavior of perturbations at a given time 
t. It is defined by 

p{i,t)^mt)\\ (26) 
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where Di(t) is a solution of equation (pT]). To determine the power spectrum, we have to 
specify the /-dependence of the initial amplitudes, Di(tin). A preferred such choice, which 



we also adopt here, is the scale invariant or Harrison-Zel'dovich spectrum [17]. The power 
spectrum is called Harrison-Zel'dovich if the variance of the mass fluctuation on horizon scales 
Rh = /o dt = t is constant, time independent: 

{i5M/M)lJ = const. 

Here (■) denotes the statistical average over many 'realisations' of perturbed universes with 
identical statistical properties. Since we know only one such universe, we assume that this 
statistical average can be replaced by a spatial average, a kind of 'ergodic hypothesis'. We want 
to express (5M)/j^(t) as a function of Di{t). Let us therefore identify the (gauge invariant) 
density variable D{x) with {6p/p){x) and let Ih be the value of / corresponding to the horizon 
size Rh = t^/Ih- Let us denote the spherical harmonics on S'^ by 3^k, where k stands for the 
multi-index m) specifying the the spherical harmonics on (see Appendix 0). We then 

3 
H 



obtain for the mass fluctuation inside a volume of size R^ 



{5M)R^{t) = / d^xh^/^5p{ic,t) = p d^xh^/\6p/p) 

JVh JVh 



P 



[ d'xh'/'j:yj^-^] (t)^p\/^EEWA(t) 

■'^H k \ P / k 1<Ih j,m 



For the final approximation, we have assumed that perturbations on scales smaller than tt/Ih 
average to zero due to the integration over Vh, and that perturbations on scales larger than 
71 /Ih are approximately constant in a volume of size Rh ~ tt/Ih, such that integration over 
Vh just gives rise to the factor Vh ( = volume of a three dimensional patch of diameter 2Rh 
on S^). With M = pVn, we then obtain 

{6M/M)l^^ Y: E 3^k(x)3^^,(x)A(t)A'W 

1,1'<Ih {j,m),{j',m') 

and 

{m/M)l^) ^ y: Di{t)Di{t) [/xh^/'y^{^)yi,{^) 

= E E E m = E ^'m)- (27) 

1<Ih j=Om=-j 1<Ih 

Let ti ^ Tc/l denote the time when the scale / crosses the horizon, Initi) = I . Since J2i<Ih ~ 
Iff, we see from (p7|) that we have to demand that 

P{l,ti = n/l) oc 

for {{5M / M)\^) to be approximately constant, i.e., for a scale invariant spectrum. 

The notion of a scale invariant power spectrum can now be used to compare the time 
evolution of density perturbations on different length scales. We want to investigate, which 
length scale collapses first. It is the scale at which the variance of the mass perturbation first 
grows of order unity. At that time, linear perturbation theory breaks down and we expect 
matter perturbations to form gravitationally bound objects. 
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3.2.2 The Final Power Spectra 



Up to an overall constant, the scale invariant spectrum is determined by the t- and l- 
dependence of the density perturbations Di{t) and by the requirement that Di{ti = tc/I) 
is proportional to l'^^"^. Then of course P{l,ti) is proportional to and the variance of the 
mass perturbation is approximately constant. 

Let us first apply this procedure to the case of radiation density fluctuations in the radiation 
dominated epoch. From solution (p!5D we find that for a fluctuation which crosses the horizon 
at times t -C vr, i.e. at times much smaller than the time of the big crunch, the maximum 
amplitude of Di is approximately constant (since the term containing cott is small) and 
therefore Ci and C2 must be proportional to /~^/^ (since a/ is proportional to to obtain the 
required behavior of Diiti). Close to the crunch, when t — vr, the expression containing 
cott diverges as (tt — t)^^. But all scales / > 2 enter the horizon at times ti < 7r/2, and 
therefore this divergence is only relevant for the mode / = 1, which enters the horizon at the 
big crunch. Of course the scales which enter the horizon already for t -C vr, also begin to grow 
as (tt — when t approaches tt. Disregarding the / = 1 mode we thus obtain 



for radiation perturbations in the radiation dominated era. At late times, close to the crunch, 
we can therefore approximate the power spectrum by P{l,t) = (?1~^{tt — This power 

spectrum obviously takes its maximum for the smallest value of /, and the induced mass 
fluctuations l^P{l) are independent of scale (see Figs. ||(A) and §(B)). 

As the next example we consider matter density perturbations in the radiation dominated 
era. We have found that for t < 1/10 they show logarithmic growth which we approximate by 
a constant. To obtain a scale invariant spectrum we therefore have to require Di oc and 
hence again P{1) oc . Only the largest scales which enter the horizon close to or after maxi- 
mum expansion do not satisfy this proportionality since we can not assume the corresponding 
density fluctuations to be approximately constant. Since the density perturbations grow with 
a certain power of t, the slope of the spectrum will decrease towards the largest scales. Matter 
density fluctuations are a special case of the composite model, when the cycle never reaches 
the matter dominated phase, a <C 1. The numerically determined power spectrum P{1) and 
the mass fluctuation l^P{T) are shown in curves A and B of Fig. |^. 

Now we determine the scale invariant power spectrum for matter density fluctuations in the 
matter dominated era. For small scales, which enter the horizon early, where Di{t) oc Ait^ scale 
invariance requires Ai oc However when the cycle approaches it's maximum expansion 
for t ^ IT, the damping term is smaller and Di grows faster, say D oc with p > 2 (around 
tmax = there is actually exponential growth, i.e. p diverges for / = 1) and we need c oc 
Towards the crunch, D is proportional to (27r — t)~^. We finally obtain roughly the following 
/-dependence of the power spectrum 




(28) 



or equivalently 




(29) 
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/>1, 0<27r-t<l 
/ ^ 1, p>2, < 27r - t < 1 
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P(l) p P(l) 




(A) I (B) 1 

Figure 4: (A): The power spectrum P{1) as a function of I for radiation density perturbations 
in the radiation dominated universe (in arbitrary units). P{1) is given at a time, when all 
scales / > 2 are inside the horizon. (B): The corresponding mass fluctuation l^P{l) as a 
function of /. 

For small / (large scales), the slope of P{1) is bigger than one (in a log-log diagram), where 
as for large I (small scales), the slope of P{1) is equal to one. This behavior is equivalent 
to the special case of the composite model with a long matter dominated epoch, a ^ 1. 
Therefore P{1) and l^P{l) in the matter dominated universe are very similar to the power 
spectra shown in Fig. ^ (E) and (F) for I < 1000. This figure actually shows a composite 
model with a = 1000. Therefore, scales with / > 1000 enter the horizon still in the radiation 
dominated era and thus do not represent this case. In a pure matter universe, there is no 
bend in the power spectrum for / ^ 1. 

Finally we approximate the power spectrum for the realistic composite model. The power 
spectrum in this case is composed of three parts. At late times, when all scales are already 
inside the horizon, we obtain (teq = T^/hq denotes the time when p^at = Prad) 

r l-\ l-^leq 

P{1) ocl I, K / < /eg (30) 

[ PP-^, 1 = 1, p>2. 

Here, the /"^-dependence is due to fluctuations which enter the horizon already during the 
radiation dominated epoch like in the fiat universe. The maximum of the power spectrum 
is expected at / ^ lf.q. In Fig. ^, (A)-(F) some examples for the power spectrum and the 
corresponding mass fluctuation are plotted. If the cycle has a long matter dominated epoch, 
the largest scale / = 1 enters the horizon soon after maximum expansion of the universe. This 
is different if the cycle does not reach the matter domination. Then the largest scale enters 
the horizon very close to the crunch and it will be the scale / = 2 which enters the horizon 
close to maximum expansion. 
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(E) ' (F) 

Figure 5: The power spectrum P{1) and the induced mass fluctuation l^P{l) as a function 
of I for matter density perturbations in the composite model (in arbitrary units). The hues 
connecting the points / G N are shown for clarity. (A) & (B): A purely radiation dominated 
cycle, a ^ 1, flmax <^ o-eq- (C) & (D): A cycle which just reaches the matter dominated 
epoch, a = 1, a^ax = ^eg- (E) & (F): A cycle including a short matter dominated epoch with 
a — A (square points, maximum of P{1) in (E) at / ~ 4) and a cycle including a long matter 
dominated epoch {amax ^ Oeg) with a — 1000 (hexagonal points, maximum of P{1) in (E) at 
I ~ 1000). 
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Figure 6: (A): The power spectrum P{1) in the composite model as a function of I, when 
matter becomes relativistic close to the crunch (arbitrary units). The cycle includes a long 
matter dominated epoch a = 1000 (The solid line simply connects the evaluated points). (B): 
The induced mass fluctuation l^P{l) for the same cycle as in (A). 

This leads to the following behavior: if a cycle is purely radiation dominated (i.e. a < aeq), 
the power spectrum takes its maximum for the smallest value of /. If a sequence of cycles 
approaches and finally enters the matter dominated era, then there will be a 'critical cycle' 
from which on the maximum of P{1) is no longer the largest scale, / = 1 but a scale comparable 

to leg > 1. 

Obviously equation (|25|) is independent of I and hence the shape of the power spectrum 
is exclusively determined by the scale invariance condition at horizon crossing and does not 
change during the subsequent growth of fluctuations. This is not the case, when the particles 
become relativistic. Then the evolution equation (^4]) does indeed depend on I and the shape 
of the power spectrum changes when the cycle approaches the crunch: further local maxima 
will occur due to the oscillating behavior of solutions of (|^), but the global maximum of the 
power spectrum remains the same. Numerical solutions for this last case are shown in Fig. |^. 

In a spatially fiat universe containing matter and radiation, the power spectrum for matter 
density perturbations can be approximated at times t ^ t^q by 



^(M)-7rTWu^- '''' 



which is similar to (|30|), only that k is continuous and the additional decrease of P for small 
values of k does not occur (In an open universe the power spectrum even starts to increase 
for the largest scales). Therefore a measurement of the power spectrum at very large scales 
(even before a cycle has reached it's maximum expansion) would (in principle) be a way to 
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1000 



Figure 7: The discrete power spectrum P{1) for the composite model with a = 1000 (square 
points connected by dots) in comparison with the continuous flat space analogon P{k) (sohd 
hne) . 



decide, whether our universe is flat, open or closed. To see the departure of (|30| ) from (^) for 
large scales, we have plotted both curves together in Fig. |^. A similar behavior has also been 
found in Wl 



3.3 Interpretation of the Results, Entropy Production in the Pre- 
vious Cycle 

We have thus found that a short time before the big crunch, the mass fluctuation, A^(Z,t) = 
l^P{l, t) is scale independent in a pure radiation universe and decreases towards large scales, 
/ < Igq in a matter/radiation universe. 

Furthermore, A^(/,t) diverges at the big crunch. Therefore, at least briefly before the 
big crunch, linear perturbation theory is no longer applicable. In the pure radiation case, we 
expect non linear effects to stop gravitational instability and prevent black hole formation 
at least on small scales. The production of gravitational entropy is thus probably not very 
signiflcant. 

However, if the universe undergoes an intermediate matter dominated period, A^(/,t) 
tends to raise towards smaller scales, approaching a very mild, logarithmic growth for / > l^g 
(see Fig. |^(D)). We also know from the corresponding flat universe analysis, if fluctuations 
grow non-linear before, due to contraction, the universe becomes radiation dominated again, 
non-linear gravity and the log-raise towards small scales will lead to the collapse of small 
scales and probably to the formation of small black holes. 

If we want to prevent excessive black hole and entropy formation in the cycle previous 
to the present one, we thus have to require that perturbations never get strongly nonlinear. 
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This yields a limit for the radiation entropy in the previous cycle. To illustrate this, let us 
assume that in the present cycle perturbations get non-linear, A^(Z,t) ^ 1 around a redshift 
of 2; ~ 10, T ~ 30K. For this not to happen in the previous cycle, we have thus to require 

T^in > 30K OtS= {\2/T^inf < 10^^- 

The radiation entropy of the previous cycle thus has to be at least a factor of 10^ times 
smaller than the present entropy. We therefore require that most of the radiation entropy of 
the universe at present, 5* > 10^'', was produced in the form of gravitational entropy from 
small fluctuations during the previous cycle. Unfortunately, we do not have a quantitative 
description for the entropy of the gravitational field (except in the extreme cases of linear 
perturbations [IS] or black holes), but it is certainly related to the dumpiness of the matter 
which is determined by the Weyl part of the curvature 0. 

We now postulate, that during the quantum gravity epoch of big crunch/big bang passage 
the entropy of the gravitational field is completely transformed into radiation entropy and the 
new cycle starts out from a state with vanishing gravitational entropy, a homogeneous and 
isotropic Friedmann Lemaitre universe. At first this postulate might seem somewhat ad hoc, 
but it is actually just what happens if a black hole evaporates due to Hawking radiation. 

It is thus feasible that most of the entropy production in the universe is actually due to mild 
gravitational clustering in the previous cycle and not due to local non-thermal processes. In 
this way it can be achieved that the radiation entropy in the previous cycle was much (several 
orders of magnitude) lower than in the present cycle. 



4 Conclusions 

We have revisited the oscillating universe and shown how it can yield a coherent solution to 
the flatness or entropy and the horizon problems of standard cosmology. We have analized 
linear gravitational perturbations in a closed universe consisting of matter and radiation. We 
can set an upper limit on the radiation entropy of the previous cycle which is at least a factor 
10^ below the entropy of the present universe. We thus postulate that most of the radiation 
entropy in the present cycle was produced as gravitational entropy by linear or mildly non- 
linear gravitational clustering in the previous cycle. During the quantum gravity big crunch/ 
big bang era, this gravitational entropy must then be transformed into radiation entropy. 

Due to the lack of of a theory of quantum gravity, we have no precise idea how this is 
accomplished. Nevertheless, this is exactly what happens if black holes evaporate! 

The reader may now object that we postulate the emergence of a Friedmann Lemaitre 
universe out of the quantum gravity era, whereas homogeneity and isotropy are naturally 
obtained in some inflationary models, e.g. chaotic inflation. However also in chaotic inflation, 
where homogeneity and isotropy is achieved by blowing up small scales, one has to require a 
cut-off of fluctuations at some very small scale, typically around Planck scale. 

We consider it to some extend a matter of taste which of the two requirements for quantum 
gravity is more 'restrictive'; that it leads to a cut-off of fluctuations around the Planck scale 
or that it leads to the transformation of gravitational entropy into radiation entropy in very 
high curvature regions. Nevertheless, it is a weakness of our model, that we cannot propose 
a clear picture of how this transformation might take place. We plan to address this problem 
in the future. 
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In our approach the creation of initial density fluctuations is not discussed. They might be 
created in (or left over from) the Planck era or they might build up during a phase transition 
in the early universe (topological defects) or by any other scale invariant process, like the self 
ordering of a global scalar field on Hubble scale. 

Also the monopole problem present in simple GUT's is not addressed. It may be solved 
by breaking electromagnetism at high temperature (see [0), or it may just not exist if the 
GUT idea of a simple unifying gauge group at high energies is not realized. 

Let us summarize the predictions and limitations of the oscillating universe in the following 
list mentioning the main problems of standard cosmology: 

• Flatness, age, entropy problem: Is solved, but predicts Q = 1 + e. 

• Horizon problem: Is solved, but a method to calculate the true age of the universe is 
still missing. 

• Cosmological constant problem: Is not addressed. 

• Monopole problem: Is not addressed, may be solved along the lines mentioned above. 

• Initial fluctuations: Are not addressed. May be left over from a quantum era or may be 
due to topological defects. 

• Fine tuning: There seems not to exist a serious fine tuning problem. But since the 
entropy was so much smaller in the previous cycle, the present cycle is the first one in 
which galaxies, stars and human beings can form. 

Clearly, the easiest way to falsify this model would be to measure f2 < 1. Many measurements 
of mass to light ratios in galaxies and clusters hint that the density parameter is indeed less 
than 1. However, these are just measurements of clustered matter. An upper Xunii on matter 
which is not or only weakly clustered on these scales is much more uncertain. Furthermore, 
if f2 = 1 we will never be able to decide whether f2 = l + eori7 = l — e, and other means to 
distinguish this scenario from, e.g., inflationary models have to be developed. 
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Appendix 



A Scalar Harmonic Functions 

In the closed universe scalar quantities like D can be expanded in a complete set of scalar 
harmonic functions 3^k(x) = 3^k(X) ^) 0) the three sphere S^: 

I^(x)=E3^k(x,^,0)^k(t), (32) 

k 

where k = (/, j, m), / = 0, 1, 2, 3, . . ., j = 0, 1, . . . , / — 1, m = — j, — j + 1, . . . , j. The variables 
X e [0,7r], 6' G [0, vr] and G [0, 27r] denote the angles describing the position on the three 
sphere. The functions J'k satisfy the Laplace-Beltrami equation with eigenvalue —k?: 

(A + e)y^ = 0. 

Here A = V-'Vj denotes the three-dimensional Laplacian on S^, k"^ = l{l + 2)K and for our 
case of interest K > (In most of the sequel we set K = 1). The harmonic functions 3^k(x) 
are given by 

3^k(x) = uff{x)Y,UO, 0), P^ = e + K={l + ifK, 

where Yjm{6, 0) are the usual spherical harmonics on and the 11^^ ■* can be expressed in 
terms of generating functions 

sin^ X ( _ 



■ i / J \ 



WAx)=^'7V7T7^{TZZZr.\ cos(/3x), 
where is the normalization factor 

n=0 

The normalization of the functions 3^k(x) is as usual 
d^'x /i'/'(x)3^k(x)X,(x) = (5k,k', 

where h{x.) is the determinant of the 3-metric of constant curvature K = 1 and (5k,k' is the 
Kronecker delta 



1, if k^(/,j,m) = k' = (/',/,m') 
0, else. 



Furthermore we choose the phases of 3^k such that 3^k(^) — 3^-k(x), with — k = {l,i^—m). 
(See e.g. or for further details). 
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B The Passage between two Cycles 



In this Appendix we want to show briefly, how the transition from a big crunch to a subsequent 
big bang can be described by an effective model. The main idea is the appearance of a signature 
change in the metric from Lorentzian to Euclidian and back. By this mechanism, the singular 
behavior of spacetime at a = disappears and the topology of the transition region is that of 



In analogy to the change of signature idea of Hartle & Hawking in quantum cosmology 



Ellis and Ellis et. al. [|T3| have shown that the classical Einstein field equations allow a 
change of signature when the metric is allowed to possess a mild singularity. The classical 
case leads to interesting possibilities for the description of an oscillating universe. 

The signature change is implemented into the metric by introduction of a lapse function 
nir): 



ds^ = -n{r)dr^ + a^{T) 



+ r\d9' + sin' 9 d^') 

1 — /cr^ 



(33) 



Here r denotes cosmic time: dr = a dt. For the discontinuous choice of the lapse function 
n{r) = e with e = ±1, there exists a surface of change S, where the metric changes its 
signature. From (^) and Einstein's equations one derives the Friedmann and Raychaudhuri 
equations for the scale factor a(r), which then hold in the regions where e = +1 and in 
V^, where e = — 1, but not on the surface of signature change S, since there the metric tensor 
is not invertible. By choosing suitable (physically motivated) jump conditions on S, one can 
find solutions for the scale factor which pass continuously through the surface of signature 
change p|, 

For the simple case of a scalar field G R with Lagrangian 



such that = (the no rolling case), Ellis et. al. have shown that the scale factor behaves for 
k=l like 

y H ^cosh(iir), r > . 

The corresponding space-time has no boundary and is geodesically complete (i.e. it has no 
singularity and the geodesies can be continued smoothly through the surface of signature 
change). Only the length of the tangent vector jumps for photons and spacelike geodesies at 
the surface of signature change E. Obviously the scale factor given by (|3^) inflates for r > 0. 
We do not have an equivalent simple example with the same nice features which does not 
inflate. However there are other (rolling) solutions to get a successful exit from inflation, but 
these solutions do not have the 'no-boundary' property of the above mentioned case. For a 
detailed discussion see []12| and |T^. Our main point is, that there is a possibility to continue 



the evolution of the universe smoothly through the crunch which might serve as an effective 
theory for the passage between two subsequent cycles. 
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Figure 8: A Geodesic of with initial velocity v^, > 0. 
B.l Geometric Representation 

It is interesting to note that only geodesies describing massive particles at rest pass through 
the point where a = 0, whereas a particle with velocity f ^ 7^ on S enters the Euclidian 
regime with an angle a as shown in Fig. ^ (The coordinates shown are t and r. The angles 6 
and (f) are suppressed, since we assume them to be constant for the indicated geodesic). 
The coordinates of the plane with angle a in the y — z plane, are 

{x,y, z) = (x, —ztana, z). 
With polar coordinates chosen as 

a; = cos ^ cos 99, ?/ = cos ^ sin (yj, z = —sai9, 
where ip G [0, tt] and 6* G [0, 7r/2 — a], one finds 

cos 6 sin if = sin 6 tan a. 
Hence a geodesic on through (f = 6 = is given by 

(f{6,a) = arcsin [tan ^ ■ tana]. 

Here 6 = Ht is the time coordinate. The velocity fs on S then is 

dip 

= 377 0=0 = tana, 
ad 

indicating that a point particle which is not at rest, does not reach the point where a = 0. A 
photon with f ^ = 1 enters with the angle a = arctanl which is 7r/4. Hence a photon is not 
represented by the boundary at a = 7r/2, but by a line inside the Euclidian regime. 
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